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Классификатор биомедицинских сигналов 


В статье рассматриваются особенности построения и принцип функционирования цифрового фильтра 
в составе нейросетевого классификатора биомедицинских сигналов. Предлагается новый подход к 
реализации процесса классификации с использованием дискриминантных функций. 


Введение 


Анализ особенностей распознавания сигналов и изображений в медицинской 
диагностике показал актуальность и перспективность дальнейшего развития и широ- 
кого использования медицинских экспертных систем на базе новейших достижений в 
области нейротехнологий и нечеткой логики [1]. При этом лечебно-диагностический 
процесс можно представить в виде процесса обмена информацией между пациентом и 
врачом, причем от пациента данные о его состоянии поступают в виде биосигналов. 
В большинстве публикаций [2], [3] приводятся примеры обработки, анализа и распо- 
знавания биоэлектрических сигналов (БЭС). Эти биомедицинские сигналы наиболее 
приемлемы для классификации из-за возможности определения большинства призна- 
ков, которые доступны в процессе наблюдения (амплитудные, частотные, фазовые), а 
также из-за наличия развитых методов и средств их обработки и анализа [2]. При этом 
специфика процесса классификации БЭС предопределяет преимущества использова- 
ния метода классификации с вычислением дискриминантных функций (ДФ) по началь- 
ным данным с ограниченными статистическими описаниями [2], [3]. 

Известно, что методы линейного дискриминантного анализа позволяют выб- 
рать проекцию пространства изображения (образа) на пространство признаков таким 
образом, чтобы минимизировать внутриклассовое и максимизировать межклассовое 
расстояние в пространстве признаков. При этом фиксируются «грубые» закономер- 
ности экспериментальных данных об объекте классификации, поскольку линейные 
модели адекватны простым геометрическим конфигурациям областей пространства 
признаков [3]. 

Целью работы является разработка и исследование нового подхода с исполь- 
зованием метода обработки векторных данных по разностным срезам для реализа- 
ции процесса классификации с вычислением ДФ. 


Классификация образов 
с вычислением дискриминантных функций 


В работе [2] приводятся обоснованные причины, по которым для задачи клас- 
сификации БЭС выбираются их частотные признаки. Известные нейросетевые клас- 
сификаторы биомедицинских сигналов и, в частности, БЭС [4], [5] представляют 
собой системы с обучением. А среди алгоритмов обучения для нейросетевых класси- 
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фикаторов распространены как известные алгоритмы (например, алгоритм обратного 
распространения ошибки с применением градиентных методов оптимизации, эври- 
стические методы обучения нейросетей, алгоритм имитации отжига [6]), так и моди- 
фицированные (например, модель обучения метричных классификаторов на основе 
вероятностно-комбинаторного подхода [7]). 

В работе [2] представлена структура классификатора, который выполняет пред- 
варительную обработку БЭС х(!) в такой последовательности: а) дискретизация по 


времени х(и) ; 6) ортогональное преобразование для выделения характерных признаков 
классификации У(); в) вычисление вектора подмножества признаков / = а. 2 } Ь 
Обучение классификатора выполняется для определения весов у; и порогов ©, 
классификации, что позволяет в дальнейшем сформировать систему ДФ =.(7), кото- 


рые определяют границы т классов образов 7, / =1.п1=1,т. Классификатор рабо- 
тает по критерию минимума среднеквадратичного расстояния вектора образа 7 и 


среднего образа 7, созданного по данным опытов и определяющего 1-й класс, кото- 
рому соответствует максимум величины определенной ДФ © (7). 


Таким образом, в процессе классификации п-мерного образа Х по т классам 
(= О | с использованием системы ДФ 2,(7) , {=1,т вида 


2 (2)=и 2+7, +...+м,2,-®,, 


И РИ (1) 


8„(2)=\2+У,22. +...+\ 2, - © 
принадлежность образа 7 классу С’, определяется по максимуму 1-й ДФ, т.е. 


тах 8: (2) = (2). (2) 


Известный цифровой фильтр в составе такого классификатора [2] реализует 
операции (1) и (2) последовательно, используя т-и блоков весовых коэффициентов 


(для формирования величины у, -2,), т блоков (п — 1) сумматоров в виде «пирами- 
ды» (для формирования ДФ =.(7)) и устройство выбора максимума (для определения 
2,(2)). При этом процедура классификации по ДФ имеет все признаки нейросетевой 


реализации классификации. Так для формирования т ДФ 2(7),...,2т(7) используется 
однослойный персептрон, который содержит т формальных нейронов с пороговой 
функцией активации (так называемые персептронные нейроны [6]), что следует из 
выражения (1). В результате один из вариантов нейросетевого классификатора по ДФ 
имеет два слоя, первый из которых представляет собой однослойный персептрон, а 
второй может быть реализован на базе сети поиска максимума с прямыми связями 
(Ееед-Еогуага МАХМЕТ) [4]. 

Итак, выбранный метод классификации БЭС по ДФ [2] в соответствии с известной 
классификацией методов распознавания образов [8] относится к методам, которые 
базируются на предсказаниях о классе решающих функций. Основным преимущест- 
вом таких методов является прозрачность математической постановки задачи распозна- 
вания, как задачи поиска экстремума. Эти методы вместе с логическими, лингвистиче- 
скими методами и методами, которые базируются на оценках плотности распределения 
значений признаков (или сходства и отличия объектов), относятся к интенсиальным 
методам распознавания [8], т.е. используют операции с признаками в отличие от 
экстенсиальных методов, которые используют операции с объектами распознавания. 
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Областью применения методов, которые базируются на предсказаниях о классе 
решающих функций, являются классы, которые хорошо разделяются, а система при- 
знаков является ортонормированной. К ограничениям этих методов можно отнести то, 
что заранее необходимо знать вид решающей функции, а в процессе распознавания 
невозможно учесть новые знания о корреляциях между признаками [8]. 


Использование обработки по разностным срезам 
при классификации образов 


В цифровом фильтре (ЦФ), который входит в состав рассмотренного классифи- 
катора [2], выполняются следующие операции: 


1) умножение каждого элемента 2, где /=1,и, входного образа 7 в виде вектор- 


ного массива данных на соответствующие коэффициенты и’, весовой матрицы УУ, полу- 
ченной в процессе обучения классификатора, с формированием векторных массивов 
(вектор-строк) взвешенных элементов вида: 


24° = {0 2... "7, Т=Ьт; (3) 

2) суммирование взвешенных элементов вида: 

а = \/;2, (4) 
в каждом векторном массиве А‘ с формированием сумм взвешенных элементов, 
которые представляют соответствующие ДФ 27) согласно выражению (1) и соот- 
ветствуют конкретным классам С7 распознавания образов; 

3) сравнение полученных сумм 2;{2) между собой для определения макси- 
мальной по величине суммы среди них, что является признаком принадлежности 
входного образа 7 в виде векторного массива к соответствующему классу Си по вы- 
ражению (2). 

В данной работе предлагается новый подход с использованием метода обра- 
ботки векторных данных по разностным срезам (РС) [9-12]. 

Итак, имеем следующие начальные данные: а) входной п-мерный образ = 
={21,...,25....2}; 6) т классов С={СИ,...,Сь...,Ст}; в) матрицу весовых коэффициен- 
тов И={и;}, где {=Ьт, 1 =Ьп; г) исходный вектор классификации Р={р;}, началь- 
ные значения элементов р; которого равны нулю. 

На первом этапе классификации по РС выполняют действия (п. 1), аналогичные 
действиям известного нейросетевого классификатора [2]. В результате имеем век- 
торные массивы в виде вектор-строк А“; (3), элементы аб; (4) которых соответствуют 
элементам конкретных ДФ (7) вида (1). Представим элементы ДФ а; как элементы 
матрицы А° (начального двумерного РС) вида 


Го 0 Г о 
т. ай А 
0__ 0 0 0 
А =| а. . а, |=|А|, (5) 
0 0 
| Ч. г-на | Ам | 


строки которой составляют элементы соответствующих ДФ (7). Дальнейшие дей- 
ствия над элементами матрицы А° (5) имеют следующий вид [10], [11]. 
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1. Определяется для каждого /-го столбца А} прямоугольной матрицы А"' раз- 
мерностью тхи, где /=|и, #=1, М, метрика — мера сходства в виде 


4 = шш А" = пита, ›Ё=т. (6) 
й 1 
В результате имеем векторную метрику вида 


ое вый (7) 


2. Выполняется коррекция элементов столбцов матрицы А"": 


=” 


аз=а,'-4,, (8) 
т.е. формируется матричный (двумерный) РС А’ = ай ь 


—— 
3. Проверяются условия для всех строк матрицы А ввиде: 


ЧА: =0, #=ЬМ, (9) 
МА, =0. (10) 


В первом случае, т.е. при выполнении условия (9), остается нулевым соответ- 
ствующий элемент р; вектора классификации Р, т.е.: 


р. =0 (0 


= 
и в дальнейшем эта нулевая строка А: из обработки исключается. Переход к п. 4. 
Во овом случае при выполнении условия (10) для последней нулевой строки 


А матрицы А“ формируется соответствующий единичный элемент р/ вектора клас- 
сификации Р, т.е. 
ре, [= Ьт. (12) 
Переход к п. 5. 


= 
4. Выполняется упорядочение матрицы А в процессе продвижения вправо к 
= — 
краю нулевых элементов ау во всех строках А; и формирование упорядоченного 


двумерного РС А‘. Переход кп. 1. 


5. Процесс завершен. 
Таким образом, последовательно будут обнулены все строки А; матрицы А’и 
соответствующие элементы р; вектора классификации Р, кроме элемента рь, который 


соответствует последней обнуленной строке А матрицы А’ в последнем Л-м цик- 
ле. Элемент р, имеет единичное значение, которое свидетельствует о принадлеж- 
ности входного образа 7, к /-му классу классификации Си=С.. 

Итак, можно сформулировать такое решающее правило классификации образов 
по ДФ: для п-мерного образа 7, который представлен в виде его взвешенных ком- 
понентов (признаков) а =" где 7; — ]-й компонент входного вектора 7, и’, вес 
вхождения ]-го компонента в 1-й класс классификации, 1=1,т; ]=.п , классифика- 
ция выполняется следующим образом: 


р, = [тах А) /=Ьт] > 2ЕС,, (13) 
где С={С,,...,С и! — множество классов; 


Р={ру,....рт} — выходной вектор классификации. 
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Таким образом, суть обработки векторных массивов данных по РС для класси- 
фикации образов состоит в использовании меры сходства между элементами а, 


векторного массива А“" в виде минимального элемента а’ ' (6) как общей части всех 
И Л 

элементов массива А’", где А’ ]-й вектор-столбец входного двумерного массива 
1 Л 


0 > 
А; М- количество этапов (циклов) обработки. Тогда разностный срез А’ вида 
1 


4 =" = ( шо" — 


р) у = 1.) 


можно рассматривать как вектор отличий векторного массива А“" в {-м цикле обра- 
1 


ботки. 

Таким образом, при классификации образов по РС вводится векторная линей- 
ная метрика сходства (7) в отличие от скалярной линейной метрики — Манхэттен- 
ского расстояния. Кстати, все известные количественные меры сходства, а именно 
расстояние Хэмминга, функция корреляции, направляющие косинусы, метрика Мин- 
ковского и все ее частичные случаи, мера сходства Танимото и другие представляют 
собой скалярные величины. То есть использованная двумерная метрика сходства (7) 
еще раз подчеркивает двумерный характер обработки по РС по всей матрице А° 
элементов ДФ в процессе классификации образов. 

В результате использование обработки по РС позволяет не только совместить 
выполнение базовых операций (1) и (2) классификатора БЭС, но и отказаться от 
формирования (накопления) ДФ 2.(7) вида (1), т.е. наоборот, уменьшать эти вели- 


чины одновременно на величину сходства т соответствующих одноименных эле- 
] > 

ментов в матричном массиве А’ в каждом цикле обработки. Такие действия по 

исключению наименьших функций ©,(7) повторяются до момента, когда остается 

последняя функция ©,(7), которая является максимальной среди всех ДФ и опре- 


деляет принадлежность образа 7 к /-му классу [10], [11]. 
В этом случае ЦФ в виде матричного процессора содержит умножитель (для 
м 1 
формирования величин у’, `2,), вычислительный блок (для формирования матриц А’) 
и блок анализа (для формирования вектора исходных сигналов Р) (рис. 1). Детальный 


сравнительный анализ аппаратных и временных затрат для предложенной структурной 
организации ЦФ (рис. 1) представлен в работе [12]. 


Вычисли- 
тельный 


блок 


Рисунок 1 — Структурная схема ЦФ 
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Анализ особенностей нейросетевой реализации 
цифрового фильтра 


Вычислительный блок и блок анализа в составе ЦФ (рис.1) можно рассматри- 
вать как карту РС, структура которой приведена на рис. 2 и которая содержит вычи- 
слительную карту в виде двумерной решетки ПЭ и карту признаков классификации в 
виде одномерной решетки ПЭ. 

Структура карты РС, приведенная на рис. 2, имеет все основания называться 
искусственной топографической картой, поскольку является моделью отображения 
признаков [13] при наличии следующих свойств: а) пространственное положение 
выходов соответствует конкретной области признаков данных, которые выделены во 
входном пространстве; 6) топология отображения содержит преобразование двумер- 
ного входного пространства в одномерное пространство признаков; в) преобразова- 
ние имеет адаптивный характер; г) в качестве критерия соответствия используется 
критерий максимума дискриминантных функций; д) в качестве меры сходства исполь- 
зуется вектор минэлементов как общих составных элементов векторных массивов; 
е) карта признаков имеет топологию одномерной решетки, которая задана выходным 
пространством. 


Рисунок 2 — Структура карты разностных срезов 


Итак, проанализируем архитектуру нейросетевого классификатора с учетом 
особенностей ЦФ в его составе по примеру классификатора с преобразованием вход- 
ного пространства [14]. ЦФ (рис. 2) в предложенном классификаторе содержит два 
слоя (не учитывая входной слой), причем количество столбцов вычислительной карты 
соответствует размерности и пространства информативных признаков. Размерность 
выходного слоя соответствует количеству классов т, к которым могут относиться 
объекты распознавания. 

У классификатора с преобразованием входного пространства [14] пространство 
входных признаков преобразовывается в пространство большей размерности, а по- 
том для классификации линейно разделимых классов в новом пространстве исполь- 
зуется простой персептрон с одним слоем обученных связей. Кроме того, между вход- 
ным и скрытым слоями находятся связи, которые не обучаются, а между скрытым и 
выходным слоями — обученные связи. Количество слов аного Ее рав- 


но трем, учитывая входной слой. атри ца элементов 
В предложенном классификаторе простран ЗНАКОВ ТО Ч $Хфун КЦИ! 


образовывается в пространство большей (двум , 


т 
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цессе обработки по РС данных вычислительной карты (скрытого слоя) через немо- 
дифицированные связи с картой признаков (выходной слой) формируется выходной 
вектор, который свидетельствует о принадлежности входного образа к определен- 
ному классу. В предложенном классификаторе, наоборот, связи между входным и 
скрытым слоями являются модифицированными (обученными), а между скрытым и 
выходным слоями — такими, что не обучаются. Учитывая наличие одной матрицы 
настраиваемых связей, как предложенный классификатор, так и известный класси- 
фикатор с преобразованием входного пространства или ВТС-классификатор [14] 
имеют определенные преимущества, связанные с обучением таких НС. 

Итак, для классификации образа 7 сформировано пространство топологических 


признаков таким образом: каждой строке А’ начальной матрицы А’ ставится в соот- 


ветствие элемент р, вектора топологических признаков Р: 
0 ие 
А; > р, =Ът (15) 


или начальной матрице А° размерностью тхи ставится в соответствие вектор 
топологических признаков Р размерностью т: 


А’ЪР, (16) 

т.е. выполняется сжатие входного двумерного пространства данных в одномерное 
пространство топологических признаков. 

При этом вектор топологических признаков Р или выходной вектор классифи- 

кации вычисляется по принципу разностно-срезовой обработки векторных данных 

(по столбцам матрицы А °) [10], [11]. Соответственно, веса нейронов изменяются по 


столбцам матрицы А"' по формуле (14), а перед обучением веса инициализируются по 
формуле (4). Используется обучение без учителя, этот метод обучения является конку- 
рентным, поскольку исключаются нейроны в строках матрицы А° с нулевой актив- 
ностью. Процесс обучения имеет итерационный характер и выполняется до момента 
стабилизации, т.е. до обнуления всех весов. Выходом НС на базе карты РС (рис. 2), 
которая фактически представляет структуру ЦФ (рис. 1), являются координаты строки 
нейронов матрицы А °, которая является победителем. 


Выводы 


1. Использование РС в процессе классификации БЭС с формированием ДФ 
обуславливает применение векторной линейной метрики сходства в отличие от 
известных скалярных метрик (метрики Минковского, расстояния Хэмминга, евкли- 
дового и Манхэттенского расстояний и т.п.). Это подчеркивает двумерный характер 
обработки по РС по всей матрице элементов ДФ в процессе классификации БЭС, что 
позволяет отказаться от формирования (накопления) ДФ с последующим их сравне- 
нием за счет уменьшения величин ДФ одновременно на величину сходства соответ- 
ствующих одноименных элементов матрицы в каждом цикле обработки. 

2. Особенность предложенной многослойной нейросети, которая используется 
в составе ЦФ для классификации образов, состоит в применении конкурентного обу- 
чения без учителя, который имеет итеративный характер с завершением в момент 
стабилизации (обнуления) строк матрицы, которая загружается в скрытый двумер- 
ный слой нейросети в начале работы. Таким образом, предложенная нейросеть имеет 
признаки карты с самоорганизацией (по аналогии с картой Кохонена), в которой 


94 «Искусственный интеллект» 32010 


Классификатор биомедицинских сигналов 2М 


выполняется сжатие двумерного пространства входных данных в одномерное про- 
странство топологических признаков, т.е. формирование вектора принадлежности 
входного образа к определенному классу. 
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